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ABSTRACT 


The primary objective of this research program is to improve our understanding of the 
dynamics and energetics of the solar corona both in the quiescent dynamic equilibrium state when 
coronal structure is dominated by the equatorial streamer belt and in the eruptive state when 
coronal plasma is ejected into the interplanetary medium. Numerical solutions of the time- 
dependent magnetohydrodynamic (MHD) equations and comparisons of the computed results with 
observations form the core of the approach to achieving this objective. Some of the specific 
topics that have been studied are: (1) quiescent coronal streamers in an atmosphere dominated 
by a dipole magnetic field at large radii, (2) the formation of coronal mass ejections (CMEs) in 
quiescent streamers due to the emergence of new magnetic flux and due to photospheric shear 
motion, (3) MHD shock fomation near the leading edge of CMEs, (4) coronal magnetic arcade 
eruption as a result of applied photospheric shear motion, and (5) the three-dimensional structure 
of CMEs. The main results from each of these studies are summarized below. 
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SUMMARY OF RESEARCH 

The formation of MHD shocks near the leading edge of CMEs is discussed in papers 1, 
2, 4, 6, 9, and 12 above. The general approach followed in this study was to begin with a 
sufficiendy simple model corona and geometry that the shock and wave formation processes 
could be studied in detail without unnecessary complicating factors. When the shocks were well 
understood in this simplified model, a more realistic corona was used. The initial studies, then, 
were carried out with a static atmosphere without gravity so the initial thermodynamic conditions 
were uniform. The use of this simple model atmosphere offered the additional advantage that 
the types of shocks that should form for various values of the parameters could be determined 
analytically. One of the main results to emerge from this study was to demonstrate that all three 
types of MHD shocks (slow, intermediate, and fast) may form near the leading edge of CMEs 
for various parametric regimes. As would be expected, slow shocks form at the lowest CME 
speeds, and fast shocks for the highest speeds. For CMEs with either slow or intermediate 
shocks, fast expansion waves travel out ahead of the shocks and produce a large enoug 
reduction in density (and therby in white-light brightness) that it should be noticeable in the 
observations. Numerical simulations confirmed the predictions of the analytic studies. Both the 
analytic and numerical work for this model predict that (1) a slow shock front should be concave 
upward (away from the solar surface), (2) a configuration containing intermediate shocks should 
also be concave upward near the center of the CME and concave downward at some distance 
from the vertical CME centerline, and (3) a fast shock configuration should be concave 

downward. 


The above studies using a static atmosphere were extended to the use of a more realistic 
coronal atmosphere containing a streamer with outflow along open field lines and no outflow 
within the closed-field region. An important feature of this streamer is that it was constructed 
so as to contain physical conditions similar to observed values, and, in particular, the plasma beta 
(ratio of thermal pressure to magnetic pressure) is less than unity throughout the streamer except 
within the narrow current sheet above the closed-field region. The streamer was computed 
numerically using a relaxation procedure developed with previous NASA support (Steinolfson 
and Hundhausen, J. Geophys. Res., 93, 14,261, 1988) and is discussed in more detail in papers 
9, 10, and 12. By examining the characteristic wave speeds in this ambient streamer, it was 
estimated that the different shock configurations predicted from the above simpler model should 
form for CME speeds within the following speed increments: slow shocks, 200-300 km/sec; 
intermediate shocks, 300-900 km/sec; fast shocks, >900 km/sec. These estimates were 
substantiated with numerical simulations of CMEs produced in the streamer as a result of an 
increase in magnetic flux within the closed-field region at the streamer base. Consequently, for 
most typical CME speeds, intermediate shocks should form near the leading edge. The 
distinctive flat-topped or concave-upward appearance of many observed CMEs may be a result 
of the formation of intermediate shocks. These more recent results are presented in an 
abbreviated form in paper 9 of the above list, and a more thorough discussion is in preparation 
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(paper 12). 

Another topic involves the simulation of the evolution of coronal magnetic structures in 
response to slow photospheric shear motion of the field line foot points, which is discussed more 
in papers 5, 6, and 11. The most important result from this portion of the research is to show 
that a single isolated magnetic arcade can be sheared slowly until it eventually erupts providing 
the coronal field lines are forced to move with the applied shear motion. In a number of studies 
by other investigators of the same problem, this condition was not applied, and the arcade did 
not erupt. In these latter studies, the field lines did not track the shear, thereby resulting in the 
build-up of what has been referred to as a resistive boundary layer at the shear boundary. For 
a typical result from the study in paper 5 the applied maximum shear speed is 2.6 km/sec or 
0.005 times the reference Alfv6n speed. In this case the sheared field moves out very little for 
the first 16 hours and then begins moving out slowly at first until it ultimately erupts outward. 
This type of behavior is precisely what is found when the outward motion of identifiable features 
in observed CMEs are followed in time. It is also shown that if the field is allowed to resistively 
slip as in other studies, the field does not erupt. 

The results in paper 5 are for a global-scale magnetic structure and demonstrate the 
feasibility of causing a single arcade to display eruptive behavior as a result of shear motion - 
providing the coronal field moves with the applied shear motion. It was speculated that the 
eruption did not occur due to an instability or a loss of equilibrium, but rather began gradually 
at first and eventually became more rapid when energy due to the shear motion was put into the 
corona (by Poynting flux through the surface) at a rate faster that the corona could adjust and 
maintain a quasi-static equilibrium. This behavior should be differentiated from a loss of 
equilibrium because there are equilibrium states consistent with the boundary conditions at each 
instant in time, but the corona simply does not have time to attain those equilibrium states before 
more energy is dumped into the corona and the boundary conditions are changed. The problem 
essentially reduces to a question of how rapidly the corona can find an equilibrium when the 
boundary conditions are changed. 

Since no further evidence (beyond a limited parametric study) was given in paper 5 to 
support the above interpretation, an attempt was made to provide some supporting results by 
more detailed examination of the shearing of active region-sized arcades similar to those used 
in previous studies by others, as opposed to the global arcades considered above. Some of the 
important results from this study are shown in the figure below (more complete results are given 
in paper 1 1), which shows the outward trajectory of a relevant field line rooted near the edge of 
the shear region for several different cases. In every case the maximum shear velocity is 0.5 
km/sec or 0.001 times the reference Alfvdn speed. For case A the shear is applied continuously, 
and, as is obvious in the figure, the field erupts outward. Although the time scales are much 
different now (with respect to the above global simulation) due to the drastically different 
physical scales (about 2 orders of magnitude smaller), the physics involved in the two studies 
should be similar. Even though the shear is removed after 10 minutes for case B, the field still 
erupts, although at a slightly slower rate. The fact that the field erupts even though the shear is 
removed suggests that a cumulative feedback between the shear and the boundary conditions does 
not artificially generate the eruption. The shear is removed after 5 minutes for case C, and the 
field does not erupt. For case D the shear is removed after 5 minutes and then reapplied after 
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brine off for 10 minutes, and the field then erupts in a manner almost identical to that for case 
A wLmrie shear is applied continuously. One intetpretation of these results is that when the 
shear applied* for a shot, enough time and the field is no. too highly stressed, the cornua t ctm 

find an equilibrium. However, when the shear is applied continuously or when ‘he firidis 
initially in a highly stressed state and the shear is then applied even for a relatively short nme, 
the corona cannot find a quasi-static equilibrium and erupts outward in an attempt to find an 
equilibrium. Note that since the rate of energy input due to the shear motion depends hne£nly 
on the sheared component of the magnetic field, even the smallest shear speed has the potential 
of eventually putting energy into a highly stressed corona so rapidly that it loses quasi-static 

equilibrium (it just takes longer). 



Additional studies discussed in paper 11 include putting in a 
aotilvins the shear in the high beta portion of the lower solar atmosphere (photosphere ). Once 
again* the shear leads to an eruption of the coronal arcade. Another study involved P ut “8 
shear viscosity term in the momentum equation as in some of the other studies of this problem. 
The viscositvdelayed the eruption, but the shear eventually drove the arcade untilit erupted. For 
problem the coronal resistivity was changed by two orders of 

Revnolds number of 10 4 to 10 2 ) at a fixed distance above the shear surface. If the field slippe 
SgL one of tire two lesistiv! regions, then this would show up as a Wnk m the fieM me^ 
the location of the changing resistivity (i.e., a resistive boundary layer would form). The 
Absolutely no indication^ die formation of a resistive boundary layer or any tank 
lines The field lines simply passed smoothly through the changing resistivity. This 
s^por, m die metiiod of applying die boundary conditions in which die field lines are forced ,o 

move with the shear motion without any slippage. 

Due to the currendy restricted observational view of CMEs, limited i to plane-of-the-sky 
images as viewed torn Earth, one of the major unresolved issues has becn ^jher tteyta™ tta 
. JL nf an e vnanding loop (arcade) with a finite extent perpendicular to the plane-of-the sky or 
ofT rapidly inflating bubble. Additionally, die change in die observed 

CME shape as it originates further from the limb is largely unknown. 
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Some insight into the 3-D structure of CMEs has been provided by the study presented 
in paper 12 for one typical case. The emergence of new magnetic flux into a finite longitudinal 
extent of a quiescent equatorial streamer belt is considered. In order to more directly relate the 
results to observations, the equations for Thomson scattering were used to convert the computed 
3-D density distribution to the white-light signatures that would be recorded m a particular 
coronagraph. The results indicate that the CME has the same basic shape in the plane of the sky 
as it occurs further from the limb, but it becomes dimmer as the distance from the limb increases. 
This simulated CME is shaped more like an expanding arcade than a bubble. 
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